EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP/2011-138 
2011/09/13 



CMS-SUS-1 1-003 



Search for Supersymmetry at the LHC in Events with Jets 
and Missing Transverse Energy 



The CMS Collaboration 



Abstract 



A search for events with jets and missing transverse energy is performed in a data 
sample of pp collisions collected at ^/s = 7TeV by the CMS experiment at the LHC. 
The analyzed data sample corresponds to an integrated luminosity of 1.14 fb . In 
this search, a kinematic variable, aj, is used as the main discriminator between events 
with genuine and misreconstructed missing transverse energy. No excess of events 
over the standard model expectation is found. Exclusion limits in the parameter space 
of the constrained minimal supersymmetric extension of the standard model are set. 
In this model, squark masses below 1.1 TeV are excluded at 95% CL. Gluino masses 
below 1.1 TeV are also ruled out at 95% CL for values of the universal scalar mass 
parameter below 500 GeV. 
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The standard model (SM) of particle physics is generally considered to be valid only at low 
energy scales and is expected to be superseded by a more complete theory at higher scales. 
Supersymmetric (SUSY) extensions to the SM [1-8] introduce a large number of new particles 
with the same quantum numbers as their SM partners, but differing by half a unit of spin. If Im- 
parity conservation |9l is assumed, supersymmetric particles, such as squarks and gluinos, are 
produced in pairs and decay to the lightest, stable supersymmetric particle (LSP). If the LSP is 
neutral and weakly interacting, a typical signature is a final state of multi-jets accompanied by 
significant missing transverse energy, J?t- Experiments at the Tevatron |[T0HT3| , SppS |T4l[T5| , 
HERA fT6llT7l[ , and LEP [18] colliders have performed extensive searches for signs of SUSY. In 
2010, the Large Hadron Collider (LHC) at CERN delivered an integrated luminosity of almost 
50 pb _1 at a centre-of-mass energy yjs = 7TeV, leading to several new searches from both 
ATLAS IHMl and CMS H23H2ZI- ' 

This Letter presents a search for SUSY based on a data sample corresponding to an integrated 
luminosity of 1.14 ± 0.05 fb _1 . The search strategy follows Ref. p3| and is designed to be sen- 
sitive to £ x signatures in events with two or more energetic jets. The search is not optimized 
for any particular model of SUSY and is applicable to other new physics scenarios with a # x 
signature. In this Letter, nevertheless, the results are interpreted in the constrained minimal su- 
persymmetric extension of the standard model (CMSSM) [ 28-30 1. The CMSSM is described by 
the following five parameters: the universal scalar and gaugino mass parameters, niQ and tni/ 2 ', 
the universal trilinear soft SUSY-breaking parameter, Aq; the ratio of the vacuum expectation 
values of the two Higgs doublets, tan /3; and the sign of the Higgs mixing parameter, p.. We con- 
sider only parameter sets for which the LSP is the lightest neutralino. The following example 
parameter set, referred to as LM6, is used to illustrate possible CMSSM yields: mo = 85 GeV, 
m 1/2 = 400 GeV, A = 0, tan f> = 10, and p > 0. 

A detailed description of the CMS apparatus can be found in Ref. |3"T fl . Its central feature is a 
superconducting solenoid providing an axial magnetic field of 3.8 T. The bore of the solenoid 
is instrumented with several particle detection systems. Charged particle trajectories are mea- 
sured by a silicon pixel and strip tracker system, with full azimuth ((p) coverage and a pseudo- 
rapidity (tj) acceptance from —2.5 to +2.5. Here, tj = — ln[tan(0/2)] and 9 is the polar angle 
with respect to the counterclockwise beam direction. A lead tungstate crystal electromagnetic 
calorimeter (ECAL) and a brass/ scintillator hadron calorimeter surround the tracking volume 
and provide coverage in tj from —3 to +3. The forward hadron calorimeter extends symmetri- 
cally the coverage by a further two units in tj. Muons are identified in gas ionization detectors 
embedded in the steel return yoke of the magnet. The CMS detector is nearly hermetic, which 
allows for momentum-balance measurements in the plane transverse to the beam axis. 

The offline event reconstruction and selection criteria described below are explained in more 
detail in Ref. [23]. Jets are reconstructed from the energy deposits in the calorimeter towers, 
clustered by the anti-Ax algorithm 1132 1 with a size parameter of 0.5. The raw jet energies mea- 
sured by the calorimeter systems are corrected to establish a uniform relative response in tj 
and a calibrated absolute response in transverse momentum pj with an associated uncertainty 
between 2% and 4%, depending on the jet r\ and pj 11331 . Jets considered in the analysis are re- 
quired to have transverse energy Ej > 50 GeV. Events are vetoed if any additional jet satisfies 
Ej > 50 GeV and \tj\ > 3, or rare, spurious signals are identified in the calorimeters |34ll35|. 
The highest- Ej jet is required to be within the central tracker acceptance and the two highest- Ex 
jets must each have Ej > 100 GeV. To suppress SM processes with genuine t x from neutrinos, 
events containing an isolated electron l|36| or muon l|37) with pj > 10 GeV are vetoed. To select 
a pure multi-jet topology, events are vetoed in which an isolated photon [38 1 with pj > 25 GeV 
is found. 
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The following two variables characterize the visible energy and missing momentum in the 
transverse plane: the scalar sum of the transverse energy Ej of jets, defined as Hj = Y^={ Ej, 
and the magnitude of the vector sum of the transverse momenta pj of jets, defined as Ht = 

| Y<i=\ Pt\/ where Nj e t is the number of jets with Ej > 50GeV. Significant hadronic activity in 
the event is ensured by requiring Hj > 275 GeV. Following these selections, the background 
from multi-jet production, a manifestation of quantum chromodynamics (QCD), is still several 
orders of magnitude larger than the typical signal expected from SUSY. While the bulk of these 
multi-jet events do not exhibit significant B j, large values can be observed due to stochastic 
fluctuations in the measurement of jet energies or mismeasurements caused by nonuniformities 
in the calibration of the calorimeters or detector inefficiencies. 

The oii kinematic variable, first introduced in Refs. Il39fj41~| , is used in the selection to efficiently 
reject events either without significant Bj or with transverse energy mismeasurements, while 
retaining a large sensitivity to new physics with genuine B j signatures. For events with two 

jets, the variable is defined as a.j = Ej ]2 / Mj = Ej ]2 / \J Hj 2 — Hj 2 , where Ej ]2 is the transverse 
energy of the less-energetic jet, and Mj is the transverse mass of the dijet system. For a per- 
fectly measured dijet event with Ej ]l = Ej ]2 and jets back to back in (p, and in the limit of large 
jet momenta compared to their masses, the value of ocj is 0.5. In the case of an imbalance in the 
measured transverse energies of back-to-back jets, olj is smaller than 0.5. Values significantly 
greater than 0.5 are observed when the two jets are not back to back and balancing genuine Bt- 
For events with three or more jets, a dijet system is formed by combining the jets in the event 
into two pseudo-jets. The total Ej of each of the two pseudo-jets is calculated as the scalar sum 
of the measured Ej of contributing jets. The combination chosen is the one that minimizes the 
Ej difference between the two pseudo-jets. This simple clustering criterion provides the best 
separation between QCD multi-jet events and events with genuine Bt- Events with multiple 
jets with Ej < 50 GeV or with severe jet energy under-measurements due to detector inefficien- 
cies can lead to values of ocj slightly above 0.5. Such events are effectively rejected by requiring 
aj > 0.55 and by applying dedicated vetoes, described further in Ref. p3). These final selec- 
tions complete the definition of the hadronic signal sample. A disjoint hadronic control sample 
consisting predominantly of QCD multi-jet events is defined by requiring ax < 0.55. 

As can be seen in Fig.[lJ the only expected remaining backgrounds with ccj > 0.55 stem from 
SM processes with genuine Bj in the final state. In the dijet case, the largest backgrounds with 
genuine Bj are the associated production of W or Z bosons with jets, followed by either the 
weak decays Z — > vv or W — > xv, where the T decays hadronically and is identified as a jet; or 
by leptonic decays that are not rejected by the dedicated electron or muon vetoes. At higher jet 
multiplicities, top quark production, followed by semileptonic weak top quark decay, becomes 
important. 

Events in the hadronic signal sample are recorded with a trigger condition that identifies can- 
didate events with energetic jets and significant Bt- Events are selected if they have Hj > 
250 GeV and Ht above a threshold that evolves with instantaneous luminosity, from 60 to 
90 GeV. In the region 275 < Hj < 325 GeV, the efficiency with which events satisfying the full 
reconstruction and selection criteria are triggered is 0.99^qq2- For events with Hj > 325 GeV, 
the efficiency is I.OO+q'qj. A set of pre-scaled Hj trigger conditions are used to record events for 
the hadronic control sample. 

The analysis makes use of two additional data samples to estimate the backgrounds with gen- 
uine Bt- First, a }i + jets sample is recorded with the hadronic trigger condition described 
above. The event selection, following closely the prescription described in Ref. ||42|, requires a 
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Figure 1: The distribution of txj, described in the text, for events in data with two or more jets 
(black dots with error bars representing the statistical uncertainties), after all event selection cri- 
teria except dj are applied and Hj > 375 GeV. For illustrative purposes only, expected yields 
from simulation are also shown for QCD multi-jet events (dotted-dashed line), associated pro- 
duction of top quarks, W, or Z with jets (long-dashed line), the sum of all aforementioned SM 
processes (solid line) and the SUSY LM6 model (dotted line). The uncertainties for the SM 
expectation, due to the limited accuracy of the available simulation datasets and jet energy 
calibrations, are represented by the hatched area. The highest bin contains the overflows. 

single, isolated muon with pj > 10 GeV in the final state and the transverse mass of the muon 
and fij system to be larger than 30 GeV to ensure a sample rich in W bosons. The muon is 
required to be separated from the closest jet in the event by Atj and Acp such that the distance 
AR = \f At] 2 + A(p 2 > 0.5. Second, a 7 + jets sample is selected using a dedicated photon 
trigger condition requiring a localized, large energy deposit in the ECAL with pj > 90 GeV 
and that satisfies loose photon identification and isolation criteria ||38|1 . The offline selection 
requires a single photon to be reconstructed with pj > 100 GeV, \tj\ < 1.45, satisfying tight iso- 
lation criteria, and with a minimum distance to any jet of AR > 1.0. For these selection criteria, 
the photon trigger condition is found to be fully efficient. 

The hadronic signal region is divided into eight bins of Hj: two bins of width 50 GeV in the 
range 275 < H T < 375 GeV, five bins of width 100 GeV in the range 375 < H T < 875 GeV, 
and a final open bin, Hj > 875 GeV. As in Ref. |23|, jet Ej thresholds are scaled down from 
their nominal values in the lowest two Hj bins to maintain jet multiplicities and thus com- 
parable event kinematics, topologies, and background composition throughout the entire Hj 
range. The background estimation methods described below are combined in the statistical in- 
terpretation of the observed data yields to provide a single prediction of the SM background in 
each Hj bin of the hadronic signal region. With respect to Ref. (23j, these refinements provide 
greater sensitivity across a broader SUSY parameter space and, in the context of CMSSM, up to 
higher-mass states. 

The p. + jets data sample provides an estimate of the contributions from top quark and W pro- 
duction (leading to W + jets final states) still remaining in each Hj bin of the hadronic signal re- 
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gion after all selection criteria are applied. Factors obtained from simulation [23J are then used 
to translate the yields in the p + jets sample to estimates in each Hj bin of the hadronic signal 
region. These factors are found to be only weakly dependent on Hj, ranging from 1.14 at low 
Hj to 0.90 at high Hj. Conservative uncertainties on all the parameters entering these transla- 
tion factors are assigned. The total systematic uncertainty is estimated to be 30%, dominated 
by the uncertainty on the efficiency for vetoing leptons. The remaining irreducible background 
of Z — > vv + jets events in the hadronic signal sample is estimated from 7 + jets events. These 
two processes have similar kinematic properties when the photon is ignored ||43l |44| , while 
the latter has a larger production cross section. Translation factors that account for the ratio of 
cross sections for 7 + jets and Z — > vv + jets, and their relative acceptances, are obtained from 
simulation ||23| and used to estimate the number of Z — > vv + jets events in each Hj bin of the 
hadronic signal region. As is the case of the ji + jets sample, these translation factors are only 
weakly dependent on Hj, ranging from 0.35 at low Hj to 0.45 at high Hj. The main systematic 
uncertainties on these factors are associated with the ratio of cross sections between 7 + jets 
and Z — > vv + jets in the simulation (30%), the efficiency for photon identification (20%), and 
the purity of the photon selection (20%), which add up in quadrature to 40%. 

Furthermore, the Hj dependence of the ratio R aj is exploited to constrain the SM background 
estimate for each Hj bin. This ratio is defined as the number of events with a.j above and 
below a threshold value of 0.55 for a given bin in Hj. The denominator of the ratio is always 
dominated by events from QCD multi-jet production and is measured in data with samples 
selected by the set of pre-scaled Hj trigger conditions. The chosen aj threshold ensures that, 
for a given bin in Hj, the numerator of the ratio is dominated by events from SM processes with 
genuine K t from neutrinos, with no significant contribution from QCD multi-jet production. 
As observed in Ref. [23], this property leads to being independent of Hj. The remaining 
backgrounds are those with genuine £j from associated production of top quarks, W, or Z 
with jets. By relaxing the aj threshold to values lower than 0.55, the numerator is instead 
dominated by mismeasured QCD multi-jet events, and an exponential dependence of R Kj on 
Hj is observed [23]. The behaviours of aj and R Kj (Hj) are observed in data and simulation to 
be robust against the effects of multiple pp collisions per beam crossing (pileup). In the statistical 
interpretation of the analysis, R Kj (Hj) is modelled as a superposition of a Hj-independent 
contribution from SM processes with genuine Ut, and an exponentially falling contribution 
to accommodate any potential QCD contamination. The latter is considered even though no 
evidence of a significant QCD contamination is found in the hadronic signal region. 

To obtain an accurate and consistent prediction of the SM background, a simultaneous binned 
likelihood fit using information from all three data samples is performed. The fit maximizes 
the likelihood L total = L hadronic x L ? , + j ets x L 7+ j ets , where L hadronic characterizes R aj (Hj) in the 
hadronic sample with a single exponential function, Ae~ kHj , to accommodate any QCD con- 
tamination and a constant, Hj-independent contribution, B, to describe SM processes with 
genuine Bj. The likelihoods L ?(+ j ets and L 7+ j ets describe the Hj-dependent yields in the yi + jets 
and 7 + jets samples. For each Hj bin, these yields are related to the numerator of R aj (Hj), 
measured in the hadronic signal sample, via the translation factors from the simulation. 

With a fit probability (p value) of 0.56, the hypothesis for the R aj dependence on Hj reproduces 
the data well, as shown in Fig. [2] The only parameter with a significant non-zero value as 
determined by the fit is the constant term B = (1.1 ± 0.2) x 10~ 5 . The other two fit parameters, 
A = (1.4 ± 1.9) x 10~ 5 and k = (5.2 ± 5.6) x 10~ 3 GeV -1 , are compatible with zero, indicating 
that no significant QCD contamination has been observed in the signal region. Furthermore, 
as a cross check, the fit is repeated with the assumption that R aj is independent of Hj, which 
in turn implies that the numerator of the ratio is fully dominated by SM backgrounds with 
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genuine B j. The result of this fit, shown in Fig.|2j has a p value of 0.41 and is in good agreement 
with the nominal fit. 
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Figure 2: The ratio R aj as a function of Hj, as measured in the hadronic data samples (black 
dots with error bars representing the statistical uncertainties). Also shown is the result of the 
simultaneous fit to the three data samples (solid line); the analogous result when assuming a 
Hj-independent hypothesis (dotted line); and, for illustrative purposes only, the expectation 
from the SUSY LM6 model when superimposed on the nominal fit result (long-dashed line). 



Table 1: Comparison of the measured yields in the different Hj bins for the hadronic, }i + 
jets and 7 + jets samples with the SM expectations and combined statistical and systematic 
uncertainties given by the simultaneous fit. 
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The fit results for all three data samples are listed in Table [TJ along with the observed yields in 
the data. Good agreement between the measured Hj distribution and the fit is observed for all 
three data samples, indicating that the observed yields are compatible with the SM background 
expectation provided by the fit. The uncertainties listed with the SM predictions are obtained 
from an ensemble of pseudo-experiments. Figure|3] compares the result of the simultaneous fit 
to the observed yields in the hadronic signal sample. 

Given the lack of an excess of events above the expected SM backgrounds, limits are set in 
the parameter space of the CMSSM. At each point of the parameter space, the SUSY particle 
spectrum is calculated with SoftSUSY [45], and the signal events are generated at leading 
order with PYTHIA 6.4 |46l . Inclusive, process-dependent, next-to-leading-order (NLO) cross 
sections, obtained with the program PROSPINO [47J, are used to calculate the observed and 
expected exclusion contours. The simulated signal events are reweighted so that the distri- 
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Figure 3: The observed event yields versus Hj in the hadronic signal sample (black dots with 
error bars representing the statistical uncertainties). Also shown are the expectations given by 
the simultaneous fit for the Z — > vv + jets process (dotted-dashed line); the associated pro- 
duction of top, W, or Z with jets (long-dashed line); the sum of QCD and all aforementioned 
processes (solid line); and, for illustrative purposes only, the SUSY LM6 model superimposed 
on the SM expectation (dotted line). 



bution of number of pileup events per beam crossing from the simulation matches that ob- 
served in data. Uncertainties on the SM background prediction, the luminosity measurement 
(4.5%) |48j, the total acceptance times efficiency of the selection for the considered signal model 
(4.5%) p3ll49| , and NLO cross section and parton distribution functions (10%) are included in 
the limit. Although signal contributions to the total yield in each of the three considered data 
samples are allowed, the only relevant signal contribution originates from the hadronic data 
sample in the case of the CMSSM. 

Figure|4] shows the observed and expected exclusion limits at 95% confidence level (CL) in the 
(mo, wii/2) plane for tan/3 = 10 and Aq = OGeV, calculated with the CL S method [50J. For 
this choice of parameters in the CMSSM, squark masses below 1.1 TeV are excluded and gluino 
masses below the same value are ruled out for tiiq < 500 GeV. The exclusion limit changes 
at most by 20 GeV in the (mo, WZ1/2) plane for different parameter values (e.g. tan [5 = 40 and 
Aq = —500 GeV), indicating that the limit is only weakly dependent on these parameters. 

In summary, the first search for supersymmetry from CMS based on an integrated luminosity 
in excess of 1 fb _1 has been reported. Final states with two or more jets and significant Bj, as 
expected from high-mass squark and gluino production and decays, have been analysed. The 
search has been performed in eight bins of the scalar sum of the transverse energy of jets, Hj, 
considering events with Hj in excess of 275 GeV. The sum of standard model backgrounds per 
Hj bin has been estimated from a simultaneous binned likelihood fit to hadronic, }i + jets, and 
7 + jets samples. The observed yields in the eight Hj bins have been found to be in agree- 
ment with the expected contributions from standard model processes. Limits on the CMSSM 
parameters have been derived and squark masses below 1.1 TeV are excluded at 95% CL in this 
model. Gluino masses in the same range are ruled out at 95% CL for mo < 500 GeV. This limit 
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Figure 4: Observed and expected 95% CL exclusion contours in the CMSSM (wiO/^i/2) plane 
(tan /3 = 10, Aq = 0, ]i > 0) using NLO signal cross sections with the CL S method. The expected 
limit is shown with its 68% CL range. The SUSY benchmark model LM6 is also shown. 

represents a tight constraint on the parameter space of SUSY models like the CMSSM. 
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